A multideterminant assessment of mean field methods for the description of electron 

transfer in the weak coupling regime 
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Multideterminant calculations have been performed on model systems to emphasize the role of 
many-body effects in the general description of charge quantization experiments. We show numeri- 
cally and derive analytically that a closed-shell ansatz, the usual ingredient of mean-field methods, 
does not properly describe the step-like electron transfer characteristic in weakly coupled systems. 
With the multideterminant results as a benchmark, we have evaluated the performance of common 
ab initio mean field techniques, such as Hartrcc Fock (HF) and Density Functional Theory (DFT) 
with local and hybrid exchange correlation functionals, with a special focus on spin-polarization 
effects. For HF and hybrid DFT, a qualitatively correct open-shell solution with distinct steps in 
the electron transfer behaviour can be obtained with a spin-unrestricted (i.e., spin-polarized) ansatz 
though this solution differs quantitatively from the multideterminant reference. We also discuss 
the relationship between the electronic eigenvalue gap and the onset of charge transfer for both HF 
and DFT and relate our findings to recently proposed practical schemes for calculating the addition 
energies in the Coulomb blockade regime for single molecule junctions from closed-shell DFT within 
the local density approximation. 

PACS numbers: 73.63.Rt, 73.20.Hb, 73.40.Gk 



I. INTRODUCTION 

Interest in electron transfer between nanoscale con- 
tacts has recently intensified, due to the advent of the 
technologically motivated field of molecular electronics 
and recent progress in experimental techniques for ma- 
nipulating and contacting individual molecules^. Elec- 
tron transport can operate between two limiting regimes, 
namely, coherent transport (CT) and Coulomb blockade 
(CB). While the conductivity is non- vanishing in the CT 
regime, CB behavior manifests itself through clear fron- 
tiers between high and low conductivity domains 3 . The 
latter was initially reported for metallic quantum dots, 
where it is dominated by a capacitive charging energy. 
On a smaller scale, CB dominated by intrinsic level quan- 
tization has also recently been demonstrated for single 
molecules between electrodes^. 

First-principle non-equilibrium Green's function 
(NEGF) methods applied to the theoretical description 
of electron transport through single molecule junctions 
are typically implemented^— in combination with den- 
sity functional theory (DFT) and a closed-shell ansatz 
(i.e., a single determinant with double occupation of the 
orbitals). This approach has been used in numerous 
works and proved to be very useful for characterizing 
the CT regime, where molecules are usually closed-shell 
systems and remain closed-shell in the electron transport 
experiments. On the other hand, the description of the 
CB regime faces additional difficulties, as demonstrated 
by Datta et al^ with the help of model Hamiltonians in 
a tight binding framework. Within this framework it was 
numerically shown that closed-shell mean-field models, 
such as DFT and spin-restricted Hartree-Fock (RHF), 
fail to yield the characteristic steps in electron transfer, 



while a spin-polarized approach such as unrestricted 
Hartree-Fock (UHF) might possibly yield a reasonable 
approximation. This is because a complete electron 
transfer generates an open-shell system containing a 
pair of singly occupied orbitals, one on the molecule and 
the other on the lead. A multideterminant configuration 
interaction (CI or Fock space ansatz) could provide a 
general solution to this problem. It is, however, difficult 
to apply a CI treatment in a first-principle description to 
the whole junction, especially since the metallic leads are 
usually represented by their mean-field band structures. 
This has been achieved at large computational costs with 
a wavefunction-based approach^, where the scattering 
boundary conditions have been formulated in terms of 
the Wigner function, making an analysis of the results 
based on molecular orbitals quite unfeasible. Within 
DFT, multideterminant schemes remain quite unortho- 
do x 12 i 13 since they diminish both its conceptual clarity 
and computational simplicity, and in addition require to 
avoid the double counting of correlation effects, which 
is far from trivial. Recently, a description of the CB 
regime was also attempted within a NEGF approach in 
conjunction with spin-polarized hybrid DF T 14 : 15 and 
Hartree Fock (HF) 16 i 17 methods for junctions containing 
finite clusters with partially filled degenerate orbitals. 

While it is evident within wavefunction theory that 
a closed-shell single-determinant ansatz is inappropriate 
for the description of an open-shell singlet state (only a 
multideterminant wavefunction can ensure spin purity), 
the situation is less clear-cut in DFT, which is in prin- 
ciple an exact ground-state theory in which a singlet 
state does not necessarily require a spin-polarized treat- 
ment. In other words, it is expected that with an ideal 
XC functional, it would be in principle possible within 
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FIG. 1: (Color online) Illustration of a molecule M exchanging electrons with the left and right reservoirs (L and R) as an open 
system where L and R are semi-infinite contacts (left panel) and where L and R are finite (right panel). 



a single-determinant Kohn-Sham (KS) DFT framework 
to treat strongly (non-dynamical) correlated systems, 
which mandatorily requires a multideterminant treat- 
ment in the framework of wavefunction theory. There 
are certainly grounds for such a belief since even stan- 
dard functionals have been shown to include implicitly a 
certain degree of non-dynamic correlatio n 18 ' 19 via their 
self-interaction effect a 20 ' 21 ; this suggests that the inclu- 
sion of a physical artefact might lead to an overall im- 
provement of the results by cancelling out the effect of 
other approximations. Following this line of thought, hy- 
brid functionals, with their decreased share of DFT ex- 
change and hence disturbed balance of error cancellation 
empirically found for local or semi-local XC functionals, 
are thus expected to perform worse with non-dynamic 
correlation issues^ 2 -. On the other hand, special func- 
tionals have been elaborated to treat bond dissociation 
within a spin-polarized KS-DFT formalism 23 . Overall, 
it is unclear how the approaches listed above, which are 
commonly used to improve DFT results, will perform in 
the context of electron transport in the CB regime. 

In this work, our goal is to gain a deeper theoretical 
insight into the requirements for a correct description of 
the charge quantization process in the CB regime by fo- 
cusing on simple model systems, which allow for a direct 
assessment of various common mean-field methods based 
on HF as well as DFT against a multideterminant refer- 
ence. To do so, we focus on a molecular triad as a model 
for an electrode-molecule-electrode system and consider 
electron transfer processes triggered by an external elec- 
tric field, either between distinct molecules or between 
two moieties connected by a bridge within a single molec- 
ular unit. Due to its simplicity, our model allows us not 
only to compare the results of mean-field versus CI ap- 
proaches on a state-of-the-art first-principle level but also 
to provide an analytical explanation for their qualitative 



differences. One of the characteristic features of the CB 
regime, namely steps in the dependence of the electron 
transfer on the applied voltage which are defined as the 
threshold voltages for charge injection, is recovered in 
our molecular model dealing exclusively with Coulomb 
effects and energy level quantization in the case of very 
weak electronic coupling between the different parts of 
the system. While a somewhat stronger coupling leads 
to a smoothening of the step-like charge transfer curves, 
we numerically show from first-principle calculations that 
the general conclusions we derive analytically for zero or- 
bital overlap still apply. 

The paper is organized as follows: in the next sec- 
tion, we introduce the general setup for our model cal- 
culations and the computational techniques employed. 
In Section [ml we derive analytically and confirm nu- 
merically why and how a single-determinant closed-shell 
ansatz fails for the case of two molecules in a cofacial ge- 
ometry with non-overlapping orbitals and why and how 
a multideterminant approach repairs the deficiencies of 
this ansatz. Section HVl is devoted to an in-depth analysis 
of the capability of standard mean-field first-principle ap- 
proaches to describe the intra-molecular electron transfer- 
process in relation to the degree of conjugation between 
the molecular donor and acceptor parts, with a special 
emphasis put on spin polarization effects. Finally, Sec.lVl 
contains a summary. 

II. METHODS AND SYSTEM SETUP 

A. Open systems vs. finite ensembles of molecules 

As a prerequisite, it is important to first argue why 
our model with three molecules should share common 
features with an open system as encountered in CB trans- 
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port experiments. The distinction between a closed and 
an open system depends on how one separates a region of 
interest from the rest of the universe 2 ^. While closed sys- 
tems obey global conservation laws for mass and energy, 
open systems in general do not. It is frequently necessary 
in the context of electronic structure calculations to par- 
tition a complex system (which might be reasonably re- 
garded as closed) into smaller components which, viewed 
individually, must be regarded as open. Transport phe- 
nomena are commonly described by differential equations 
for finite objects, with the "openness" of the system de- 
fined by the boundary conditions applied to these equa- 
tions. This can be for instance achieved by employing 
periodic boundary conditions, which are adapted to the 
requirements of linear-response theory^ or by defining a 
system which is coupled to two or more ideal reservoirs 
of particles^. In the latter case, the conductance is then 
expressed in terms of the quantum-mechanical transmis- 
sion coefficients of the system in between the reservoirs, 
which is the basis of the non-equilibrium Green's function 
formalism (NEGF) 27 for DFT- based electron transport 
calculations^. 

Our definition of "openness" is more elementary in 
the sense that the central molecule in an ensemble of 
three has a variable number of electrons (mass is not 
conserved). By applying an external field, the electrons 
can be transferred to the neighboring molecules, which 
are thus playing the role of reservoirs (Fig. [I]). This is 
related in spirit to the statistical interpretation of frac- 
tional occupation numbers in DFT, whose physical mean- 
ing is to describe time averages of electrons exchanged 
between two open systems^. The simplicity of our ap- 
proach does not allow for the definition of a conductance 
or current, or for the characterisation of screening ef- 
fects ^ 29 ' 30 . The reasons for that are two-fold: i) The 
absence of periodic boundary conditions for our "elec- 
trodes" prevents the unambiguous definition of a density 
of states or Fermi energy, which are necessary for deriv- 
ing a conductance within the Landauer-Biittiker formal- 
ism and for a steady-state description of the current; ii) 
in order to capture the electrostatic phenomenon leading 
toscreening correctly, the metallic character of electrodes 
(not present in our model) would be essential. 

Due to these limitations of our model, we cannot re- 
ally claim that we describe a CB-type transport sce- 
nario, where the key experimental quantity would be cur- 
rent/voltage (I/V) curves, which we are unable to pro- 
duce explicitly. There are, however, spectroscopic capac- 
itance experiments 3 ^, which equally demonstrate gate- 
induced charge-quantized electron transfer without re- 
lying on a direct evaluation of I/V measurements but 
instead by observing discrete peaks in the device capac- 
itance of a very sensitive transistor. Such experiments 
are much closer related to our model than the direct mea- 
surement of electron transport; their relation to Coulomb 
blockade phenomena has been clarified by Biittiker et 
al. 32 A common feature of both experimental approaches 
is the occurrence of either distinct steps in I/V curves or 



equivalently sharp peaks in its first derivative which are 
also found in the gate-dependence of the capacitance in 
Ref.— . We show in our work that a correct description 
of such steps can be achieved by a multideterminant ap- 
proach and use this solution as a benchmark for assessing 
the reliability of common meanfield techniques. 



B. Computational approaches 

In this work, we have used several computational tools 
as a basis for our argumentation. The Austin model 1 
(AMl)^ 3 -, as implemented in AMPAC^i, is a wavefunc- 
tion based parameterized technique, which makes it com- 
putationally very efficient for systems for which reliable 
atomic parameters are available, i.e., for the atoms typ- 
ically found in organic molecules but only for a rather 
small selection of metallic atoms. The conceptual sim- 
plicity of this method and its minimal basis set make 
it complementary to the analytical wavefunction based 
models developed in Sec. Mil In addition to closed- 
shell calculations, AMI has also been used in the frame- 
work of a complete active space configuration interac- 
tion (CASCI) scheme, in which all excited configurations 
among the specified occupied and vacant RHF MOs are 
included in the wavefunction expansion, and their par- 
ticipation is determined variationally. 

In contrast, the ab initio techniques used within the 
GAUSSIAN package 3 ^ do not require parameters fitted 
to experiments, which has the advantage that the issue 
of transferability of such parameters to different systems 
or different physical boundary conditions never arises; 
however, this comes at the cost of a much higher com- 
putational effort. This code allows for a direct compari- 
son (i.e., within the same computational setup and using 
the same basis sets) of density functional theory-based 
techniques, where local/semilocal or hybrid functionals 
are used for the exchange correlation part of the Hamil- 
tonian, with single and multideterminant implementa- 
tions of wavefunction theory, which is the main focus 
in Sec. IIVI The multideterminant ab initio wavefunc- 
tion calculations are performed within a complete active 
space self-consistent field (CASSCF) scheme, in which - 
as a contrast to CASCI - the linear coefficients for the 
expansion of both the wavefunction into Slater determi- 
nants and the molecular orbitals into atomic orbitals are 
optimized simultaneously. Both CASCI and CASSCF 
active spaces include 2 highest occupied and 2 lowest un- 
occupied MOs. 



C. Model junction with three molecules 

In Fig. [H we present the results of calculations per- 
formed at both the AMI /RHF and AM1/CASCI levels 
for the field-induced charge transfer in a TCNQ-TTF- 
TCNQ cofacial stack (with TTF - tetrathiafulvalene - 
a strong electron donor and TCNQ - tetracyanoquin- 
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created by so-called sparkle charges— that are shown to 
be equivalent to a homogeneous field in the next section. 
The geometry of each molecule was relaxed individually, 
whereas the impact of inter-molecular interactions and 
polarization due to sparkle charges or external field on 
the positions of atomic nuclei have been neglected, be- 
cause the focus of our article is on electronic effects only. 
The molecule in the middle (M) has been chosen as a 
donor and the molecules on the left (L) and right (R) 
sides as acceptors in order to reach a scenario where a 
significantly larger field is needed to induce the second 
charge transfer (M— >R) compared to the first (L— >M). 
In the CI description (Fig. [2}d), a series of full-electron 
transfer steps is found while the process is continuous 
with RHF (Fig. [2fc) up to the point where two electrons 
have been exchanged between L and R; in the latter case, 
only a fractional charge is localized intermediately on 
M. The unphysical absence of stepwise electron trans- 
fer in RHF is due to the inherent double occupation of 
the MOs in this theoretical framework regardless of the 
strength of the external field. Fig. [2] also shows that the 
very onset of continuous charge transfer within RHF co- 
incides with the field strength required to induce a full 
one-electron transfer within CI. We will rationalize this 
deep correspondence by means of the analytical model 
introduced in the next section. If the results in Figs. [2] 
b and c were taken from experimental current /voltage 
(I/V) curves instead of theoretical calculations on elec- 
tron transfer in an ensemble of molecules, one would in- 
terpret one-electron steps (as found with CI) as an indi- 
cator for the CB regime and would associate the mono- 
tonic evolution (as obtained from RHF) with CT. We 
thus reach the conclusion that the indiscriminate appli- 
cation of the closed-shell ansatz beyond its applicability 
may lead to a physically wrong assignment of the trans- 
port regime, which motivates a deeper analysis of the 
differences between the RHF and CI results. 



III. ELECTRON TRANSFER BETWEEN TWO 
NON-OVERLAPPING MOLECULES IN 
WAVEFUNCTION THEORY 



FIG. 2: (Color online) Intermolecular charge transfer in 
a frozen cofacial TCNQ-TTF-TCNQ stack (a) induced by 
"sparkle" charges of opposite sign (±qSp) as obtained 
from (b) Complete Active Space Configuration Interaction 
(CASCI, with the active space including the highest 3 and 
lowest 3 MOs) and (c) RHF calculations with the AMI Hamil- 
tonian. TCNQ = tetracyanoquinodimethane (acceptor), TTF 
= tetrathiafulvalene (donor). 



odimethane - a strong electron acceptor). Although sim- 
ilar results can be obtained with first-principle methods, 
we generate them here with the parameterized Hartrce- 
Fock AMI Hamiltonian^i, since the use of a minimal ba- 
sis set facilitates our analysis in the following section. 
For all AMI calculations, the external electric field is 



A. A closed-shell single-determinant approach 

When we remove the central molecule M and com- 
pute the electron transfer directly from L to R, we find 
the same qualitative differences between CI and RHF 
(Fig. [3^) ■ We therefore concentrate hereafter on this sim- 
pler system in order to contrast the results obtained by 
quantum-chemical calculations within the RHF approach 
and within the CI formalism and to rationalize the dis- 
crepancy with the help of an analytical model. 

In the field-induced charge-transfer process within the 
TCNQ-TCNQ dyad (Figs. and b), one electron is 
transferred from the highest occupied MO (HOMO) of 
the molecule playing the role of the donor, to the 
lowest unoccupied MO (LUMO) of the molecule acting 
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FIG. 3: (Color online) a) Charge transfer in the TCNQ- 
TCNQ cofacial dyad (the molecular geometry is shown in 
the inset ) induced by sparkles as obtained from AM1/RHF 
(open symbols) and AM1/CASCI (full symbols) calculations;, 
b) general energy diagram obtained via a two-state model for 
non-overlapping molecules, see text for details. 



as the acceptor, \a). In this one-electron picture, it is 
sufficient to consider only the HOMO \h) and LUMO \l) 
of the total system and the field-evolution of the weight 
of the \d) and \a) levels since the nature of the other or- 
bitals of the complex does not vary with the field. The 
calculations show that the field-evolution of the RHF 
eigenenergies for \h) and \l) can be divided into three lin- 
ear regions I, II, and III (see Fig. [3Jd). In region I, there 
is no charge transfer yet and the two orbitals \h) and 
|Z) correspond to the pure HOMO and LUMO levels of 
the individual molecules: \h)—\d), \l) = \a). In region III, 
the two-electron charge transfer is complete, and |ft,) = |a), 
|Z) = |d). In the intermediate region II, the gradual charge 



transfer found within RHF is due to a continuous mix- 
ing of the individual \d) and \a) orbitals in the frontier 
MOs \h) and |Z) of the dyad. We discuss below the ori- 
gin of this gradual orbital mixing, which in spite of the 
strictly zero orbital overlap between \d) and \a) ensured 
by the considerably large distance between the molecules, 
is characterized by a clear onset at the RHF level when 
moving from region I to II. Note that this mixing yields 
an energy separation between the \h) and |Z) levels of 
about 1-2 eV throughout the entire region II of Fig. [3K- 
An explanation could be sought in terms of an avoided 
crossing scenario, as it was done for the formally similar 
case of molecular dissociation^; however, this explana- 
tion was disproved for the latter case^ on the grounds 
that avoided crossings would require the size of the in- 
teraction matrix element to match the energy splitting, 
which in general is not the case for large intermolecular 
separations. 

In order to rationalize these results, we consider a sim- 
ple model encompassing at zero field a doubly occupied 
level on a donor D and a non-overlapping unoccupied 
level on an acceptor A, which we refer to as a Z22-modcl 
(standing for zero-overlap, 2-electron and 2-orbital). In 
this framework, the individual RHF eigenenergies of the 
isolated molecules are defined by: 



(d\ H D \d) 

(a\ H A \a 



(dd\dd) 



= (d\T + V D \d) 
= (a\T + V A la) 



{dd\dd) 
(1) 



where T is the kinetic energy operator, Vrj and V a are 
the potential energy operators describing the interaction 
with the other electrons and with the nuclei on the re- 
spective molecules; (dd\dd) is the two-electron Coulomb 
integral that appears only for the occupied orbital. When 
D and A form a system with zero overlap between the 
fragment orbitals \d) and |a), the composition of the or- 
thonormal HOMO and LUMO orbitals \h) and \l) of the 
dimer are 



\h) = cosQ \d) + sinQ \a) 
\l) = -sinO \d) + cosQ \a) 



(2) 



with the mixing parameter O characterizing the system 
completely. The RHF total energy defined only by the 
doubly occupied level \h) then becomes^: 



E RHF = 2 (h\ T + V D + V A + V bias \h) + (hh\hh) = 
2e° d - (dd\dd) + 2qn D + 2sm 2 9(e° -e° d - {dd\aa) 
-q(KA + Kn)) + sm 4 0((<ie?|dd) + (aa\aa) — 2(dd\aa))(3) 



The external field enters Erhf via Vbias which depends 
on the magnitude of the sparkle charges q (see Fig. [2] for 
details) and on the spatial distribution of the occupied 
orbital via an effective coefficient k, that is V{,i as (A,D) 

There is no mixing of the D and A fragment orbitals 
if sin 2 0=O, i.e., \h) = \d) and \l) = \a) (region I in Fig.[3b). 
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Therefore, the degree of RHF mixing can be obtained by 
minimizing ^irhf m Eq. [5] with respect to sin 2 ©, where 



i 2 = 



(dd\aa) - (e° a - e°) + (n D + n A )q 
((dd\dd) — (dd\aa)) + ((aa\aa) — (dd\aa)) 



(4) 



The condition sin 2 >0, which entirely determines the 
onset of RHF mixing, is satisfied only when the numera- 
tor of Eq. [5] is non- negative (the denominator is positive 
as one-center Coulomb integrals in general exceed two- 
center integrals), that is when 



(e° - n A q) - (erf + KDq) < {dd\aa). 



(5) 



Eq. [5] means that the onset of mixing (the point at the 
crossing of region I and region II in Fig. [3]b) takes place 
when the energy difference between the LUMO of the ac- 
ceptor and the HOMO of the donor in its linear decrease 
under the influence of the sparkle charges q, equals the 
Coulomb integral (dd\aa). Within RHF the cigcncncrgies 
of the vacant and occupied MOs approximate the elec- 
tron affinity EA and ionization potential IP, respectively, 
due to Koopmans theorem. On the other hand, the inte- 
gral (dd\aa) is the Coulomb interaction energy between 
the charge densities of the two molecules, which are well 
separated in space in our case. This Coulomb integral is 
expected to evolve as 1/d for the asymptotic limit of large 
intermolecular distances, which is verified in Fig. U show- 
ing the distance dependence of the HOMO-LUMO gap of 
the dyad at the onset of the charge transfer as calculated 
numerically with AMI. Therefore, in the long distance 
limit, the RHF mixing onset is also the point where the 
difference between EA and IP is fully compensated by 
the Coulomb attraction of the formed ions, thereby gen- 
erating the necessary conditions for a full one-electron 
transfer, which occurs within the correct CI treatment 
at the same bias in Fig. [3k.. 

It is straightforward to demonstrate from the standard 
expression of the RHF energies of the occupied and va- 
cant molecular orbitals^: 



e h = (h\ t + V D + V A + V bl as | ft) + (hh\hh) 
€l = (l\t + V D + V A + Vbias 10 + 2{hh\ll) + [hl\hl) (6) 

and making use of Eqs.[2]and[I]that the HOMO-LUMO 
gap actually remains constant and equal to (dd\aa) 
throughout the mixing region II, though the eigenener- 
gies eh and ei may change. 

We note further that the Mulliken charge^ Q on the 
non-overlapping fragments D and A is proportional to 
the square of the coefficients in the orbital expansion of 
\h) multiplied by |ft.)'s occupation (which is always 2), 
that is |Q|=2 sin 2 ©. According to Eq. HJ Q is linear 
with the external field generated by the sparkle charges 
q, which explains the continuous intermolecular charge 
transfer obtained from closed-shell calculations (Fig. [3k). 
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FIG. 4: (Color online) Evolution of the threshold 
HOMO/LUMO gap for the onset of charge transfer in the 
AM1/RHF calculations (i.e., region II in Fig.[3b) in a TCNQ- 
TCNQ cofacial dyad as a function of the intermolecular 
distance d. The diagonal line shows the point-charge in- 
teraction energy, evolving as E=f/d (in atomic units e=f, 
lHartree(Ha) = 1 /Bohr) . 



B. Multideterminant approach 

The open-shell singlet wavefunction describing the fi- 
nal state of an electron transfer from D to A can be writ- 
ten as 

f 



*(D+A") = -^(|da| - M) 
= -l=(da + ad)(U - Tl) 



(7) 



where \da\ is a Slater determinant, and (da + ad) = 
(|d(l)> |o(2)> + |o(l))|d(2)» and (U - IT) = 
(|t (1)) |I (2)) - || (1)) |t (2))) arc the shorthand nota- 
tions for the spatial and spin parts of the two-electron 
wavefunction, respectively. In the following, we illustrate 
how to construct from the RHF-MOs \h) and |Z) intro- 
duced in the last section or in other words how to move 
from a single determinant (RHF) to a multideterminant 
(CI) description, where the Fock space is a minimal CI 
space encompassing four determinants generated by sin- 
gle and double excitations from \h) to \l). 

It is straightforward to show that, in the Fock space 
spanned by four Slater determinants, namely, \hh\, \ll\ 7 
\hl\ and \lh\, a multideterminant singlet wavefunction of 



7 



0.8 




2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 

Sparkle charge q Sp , |e| 

FIG. 5: (Color online) CI coefficients x and y appearing in 
Eq.[8]as obtained from AM1/CASCI calculations for a TCNQ- 
TCNQ cofacial dyad as a function of the sparkle charge in the 
RHF mixing region II (see Fig. [3}- Inset: y vs. x as obtained 
from these calculations compared to the prediction of the Z22 
model. 



polarity, whereas the energy of D a A° is unaffected by q 
since the out-of-plane polarizability of the planar TCNQ 
molecule is negligible. The crossing of the ground and 
first excited states is not avoided in the zero wavefunc- 
tion overlap limit. We will show in the next Section that 
the abruptness of this transition is smoothened when in- 
troducing the overlap between the donor and acceptor 
orbitals. 

It is clear from Eq.[5]that the CI coefficients x and y are 
determined by the external field parameter q via the or- 
bital mixing parameter 0. In order to obtain the details 
of this dependency, we note that cos(2<d) = 1 — sin 2 Q 
decreases linearly with the external field (according to 
Eq. |4l sin 2 Q is linear with q in the RHF orbital mixing 
region II). Therefore, it follows from Eq.[5]that the coeffi- 
cient ±x varies linearly with q between -l/v2 and l/v2 
and its counterpart ±y is defined by y = \J\j2 — x 2 . In 
Fig. [5] we plot the coefficients x and y as a function of 
the parameter q for the external field, as calculated at 
the corresponding AM1/CASCI level for the TCNQ co- 
facial dimer within the mixing region found in Fig. 3a. 
This data validates the functional dependencies derived 
from our analytical model. 



the general form 



Discussion 



$(D + A-) = x(\hl\ - \lh\) + y{\hh\ ~ \ll\) 
= [(-xsin(2Q) + ycos(2Q))(dd - aa) 
+ (xcos(2Q) + ysin(2<d)){da + ad)](U - U) 



(8) 



with the MOs \h) and \l) expressed as in Eq. [21 corre- 
sponds to "F in Eq. for the whole range of the mixing 
parameter 0, if the conditions 



x = —j=cos(2Q) , y — —j=sin(2Q) 



(9) 



are fulfilled. We stress that ^5>(D + A~), in order to remain 
constant according to Eq. has to vary continuously 
when expressed in the basis of the molecular orbitals \h) 
and |Z), which are variables themselves due to RHF mix- 
ing. For the mixing parameter at the onset and in the 
middle of the region II of Fig. [31 we find by combining 
Eqs.[5]and M 



<$(D+A-) = (l/v/2)(|«|-|ift|) 
m{D+A-) = (l/V2)(\hh\-\ll\) 



for -> 

for© = 7r/4.(10) 



Whereas the RHF determinant \hh\ dominates the 
wavefunction before the RHF mixing (region I), the CI 
wavefunction does not contain a relevant contribution 
from it just at the onset (see Eq. LT0|) . This abrupt change 
in the nature of the electronic ground state, pointing to 
the inadequacy of RHF to provide a physically correct so- 
lution to the problem, is due to the crossing of D°A° and 
D + A~ at the threshold value of q. It is mainly the total 
energy of D + A~ that strongly depends on q due to its 



Our analytical model demonstrates that a closed-shell 
single-determinant RHF approach wrongly predicts con- 
tinuous charge transfer between non-overlapping moi- 
eties, which is automatically corrected to integer elec- 
tron transfers with a minimal multideterminant ansatz. 
While RHF finds variationally the best wavefunction in 
the form of a single determinant based on doubly occu- 
pied MOs, this solution is totally inappropriate for the 
description of the open-shell singlet state reached by sin- 
gle electron transfer processes. Such an open-shell singlet 
state can only be described by a series of determinants 
based on RHF MOs. The coefficients in the expansion in- 
troduced in Eq. [5] add the additional variational degrees 
of freedom to the model, so that the resulting multideter- 
minant ansatz properly describes the system. Our Z22 
model explains: (i) the sharp onset of continuous charge 
transfer in RHF as a function of the external field (which 
is different from a conventional avoided crossing) ; (ii) the 
reason for the unphysical linearity of the field-dependence 
of the charges on the donor and acceptor units in RHF; 
(iii) the origin of the HOMO-LUMO gap in the RHF 
mixing region; (iv) why the RHF charge-transfer onset 
coincides with the full one-electron transfer in CI; and 
finally (v) how CI corrects for the described artefacts of 
RHF. 

As the energy splitting between the two frontier MOs 
\d) and \a) gets reduced with the increase in the elec- 
tric field, the system as a whole becomes strongly corre- 
lated, or in other words non-dynamic correlation becomes 
essentia l 19 i 40 ' 41 ; the latter is also frequently referred to 
as static, left-right or near-degenerate correlation, and 



is also crucial for a correct description of the transition 
point in molecular dissociation processe s 42 ! 43 . There are 
similar methodological problems in the case of the dis- 
sociation of ionic molecules, e.g. LiF, where deficiencies 
in RHF or DFT calculations lead to (partially) charged 
atoms at large separation, while a correct treatment does 
always yield neutral open-shell atoms after the crossover 
in the total energies of the neutral and ionic state o 44 i 45 . 

Our results were derived so far for the case of zero 
orbital overlap between the two molecules. In the next 
section, we allow for a weak overlap between the orbitals 
of the donor and acceptor moieties and show that this 
does not alter the general conclusions derived here. 



IV. MEAN FIELD APPROXIMATIONS FOR 
INTRAMOLECULAR ELECTRON TRANSFER 

A. Aviram Ratner molecules 

After having compared single- versus multideterminant 
approaches for a simple system both at an analytical 
and semi-empirical level, we devote the current section 
to approaching the same problem with common mean- 
field techniques based on a first-principle description; all 
calculations reported hereafter were performed with the 
Gaussian package 3 ^. We have investigated different model 
systems to vary the coupling strength between the donor 
and acceptor moieties determined by the degree of conju- 
gation provided by the bridging unit and assess its influ- 
ence on the nature of the electron transfer (i.e., whether 
it is continuous, stepwise or intermediate in its depen- 
dence on the external electric field is not straightfor- 
ward to predict a priori). In particular, we consider the 
Aviram-Ratner molecule—, where the acceptor TCNQ 
and the donor TTF units are linked by a non-conjugated 
bridge denoted TCNQ-PHP-TTF, and its shorter analog 
with a conjugated bridge denoted TCNQ-PP-TTF. For 
these systems, the applied field is not perpendicular to 
the molecular planes of the D and A moieties as in the 
previous section but on the contrary directly along the 
molecular backbone so that the gradual polarization of 
the 7r electron system is also expected to play a role in 
the physics of the electron-transfer process. 

In the absence of any external field, the frontier MOs 
are completely localized on the D and A parts for the 
longer (PHP) bridge (Fig. E}. For the shorter (PP) 
bridge, there is some spill-over in the HOMO and LUMO 
orbitals with RHF and even appreciable derealization for 
DFT using semi-local XC functionals (hybrid functional 
yield a behaviour intermediate between semi-local DFT 
and RHF). As expected, we also note that the HOMO- 
LUMO gap is much larger in the RHF than in the DFT 
description providing a negligible value. Within standard 
DFT, we are dealing with a peculiar artefact for this sys- 
tem in the sense that the strongly reduced gap of the 
isolated TTF and TCNQ molecules makes the LUMO 
of the acceptor lie lower than the HOMO of the donor. 
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a) 




FIG. 6: (Color online) Shape of the frontier orbitals (HOMO 
in the left and LUMO in the right panels, respectively) at 
zero field for TCNQ-PP-TTF as obtained with a) RHF; b) 
DFT-BLYP and c) for TCNQ-PHP-TTF with RHF using in 
all cases a 6-31G** basis set. The size and color of the balls 
are representative of the amplitude and sign of the LCAO 
coefficients, respectively. 



This can only be remedied by hybrid functionals with a 
sufficient share of HF exchange (e.g. 50% in the hybrid 
BHandHLYP potential) 4 !. 

We stress that the intra-molecular electron transfer 
falls in the weak coupling or Coulomb blockade regime for 
the two molecules studied, although they differ by their 
internal degree of coupling provided by the PP and PHP 
bridges. Low thresholds and smooth charging curves 
are physically correct observables in the strong coupling 
regime of coherent electron transport since CASSCF 
calculations predict there a continuous partial electron 
transfer reflecting the change in the spatial distributions 
of the strongly overlapping (and hence delocalized) fron- 
tier orbitals. This scenario is, however, beyond the scope 
of this paper since our model derived in Sec. Mil is based 
on the assumption of zero orbital overlap, which makes 
it unsuitable for predictions in the coherent transport 
regime. 



B. Single-determinant closed-shell methods: 
non-spin-polarized DFT and RHF 

Fig. shows the charge separation between the donor 
and acceptor parts in the TCNQ-PP-TTF molecule as 
obtained from first-principle calculations within DFT at 
the local density approximation (LDA) level and with 
the hybrid exchange correlation (XC) functional B3LYP 
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FIG. 7: (Color online) Intramolecular charge transfer between 
the donor and acceptor moieties in TCNQ-PP-TTF induced 
by sparkle charges (dashed lines) or an external electric field 
(dotted lines) as obtained by ab initio closed-shell pure DFT 
at the LDA level (circles), hybrid DFT at the B3LYP level 
(squares), and RHF calculations (triangles) with a 6-31G** 
basis set. The significant zero- field charge transfer between 
the moieties obtained in the LDA calculations originates from 
an artefact: for the isolated molecules, the LUMO of the ac- 
ceptor TCNQ lies much lower than the HOMO of the donor 
TTF (-5.44 versus -3.90 eV) unlike RHF results (-1.97 and 
-6.77 eV); note that the gas-phase EA of TCNQ and IP of 
TTF, 2.8 annd 6.8 eV, respectively**. 



(with 20% of HF exchange) as well as with RHF for a 
direct comparison with wavefunction theory. The data 
in Fig. [7] further demonstrates the complete equivalence 
of the polarization induced by sparkle charges and by the 
application of an external electric field, thus justifying the 
use of sparkle charges throughout this work. We display 
in Figs. [S] and O results calculated from wavefunction the- 
ory (Figs. [SJi, ISM and DFT (Figs. [3b, [5b) for a molecule 
introducing a strong coupling between the D and A units 
(Fig. H TCNQ-PP-TTF) and a weaker coupling (Fig. H 
TCNQ-PHP-TTF) , respectively. In these figures, we also 
contrast the closed-shell results with calculations allow- 
ing for spin polarization (see next section) at both the 
RHF and DFT levels and a coupling with a CI scheme 
in the case of RHF (Figs. [5^, [5^l). The general features 
of the ab initio closed-shell RHF curves for intramolecu- 
lar electron transfer in Figs.[7l[8^ and[9^i are reminiscent 
to those encountered in the model dyad with AMI (see 
Fig. [3J . Whereas there is some smoothening at the onsets 
due to orbital overlap in TCNQ-PP-TTF (Figs. |8Ji), 



the absence of such overlap in TCNQ-PHP-TTF (Fig.Ui) 
makes the mixing region as sharply separated from the 
integer-charge regions as in the model systems. A step- 
like behavior is found with CI and a continuous change of 
the Mulliken charges with closed-shell RHF, as observed 
in Figs. [3 and H The closed-shell DFT-LDA charging 
curves are fully linear with the field, which is due to 
the possibility of fractional charges in the framework of 
DFT; using the hybrid B3LYP or BHandHLYP function- 
als with 20% and 50% of HF exchange, respectively, does 
not modify this behavior. 

It has been recently suggested that a standard DFT 
framework has inherent problems for describing electron 
transport in both the CB and CT regimes**, due to the 
self-interaction of electrons** in a Kohn-Sham framework 
and the lack of a derivative discontinuity (DD)** in the 
evolution of the HOMO eigenenergy as a function of its 
occupancy, which can be fractional. In Ref4*, these two 
issues have been portrayed as intimately linked and a SI 
correction scheme has been devised as a remedy. In par- 
ticular, Ref4* shows two distinct things: i) calculations 
with a model discontinuous potential lead to the recovery 
of steps in the electron occupation in a model quantum 
dot in the weak coupling regime, while a continuous po- 
tential induces a continuous dependence on voltage, thus 
demonstrating the role of DD; ii) effective tight-binding 
transport calculations on molecular system contacted by 
two probes suggest that the corrections for SI could serve 
to introduce the desired discontinuities into the XC func- 
tional and thereby produce CB steps. In a more recent 
paper from the same group**, it is argued that many fail- 
ures of standard DFT can be traced back to SI, though 
the authors also concede that the physical reasons for the 
lack of DD in common XC functionals are too complex 
to be explained by SI alone. It is therefore an open issue 
to assess how much of the problems DFT faces in de- 
scribing transport in the CB regime can be attributed to 
SI; note that Datta** has also stressed that SI is equally 
absent in RHF and UHF, therefore not offering an expla- 
nation why the latter performs better than the former. 
We showed above that the incorrect description of steps 
in the charge/voltage curves related to a lack of DD in 
Ref4* is actually a general feature of closed-shell meth- 
ods. Since RHF is by definition completely SI free, the 
RHF evolution in Figs. [3j [3 [8^, and [9^ exemplifies the 
limits for the improvements that can be achieved by SI 
correction techniques within non spin-polarized KS-DFT, 
at least with commonly used XC functionals. 



C. Single-determinant open-shell methods: 
spin-polarized DFT and UHF 

Spin-unrestricted or spin-polarized calculations have 
been performed in the past for a wide variety of open- 
shell systems in quantum chemistry. Within wavefunc- 
tion theory, it is well known for homolytic bond dis- 
sociation that although a multireference ansatz is the 
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FIG. 8: (Color online) Intramolecular charge transfer between 
the donor and acceptor moieties in TCNQ-PP-TTF induced 
by sparkle charges as obtained by: a) ab initio wavefunc- 
tions; and b) hybrid DFT with a 6-31G** basis set. In 
a) restricted (circles), PO-unrestricted (triangles) and BS- 
unrestricted (squares) HF and CASSCF (diamonds, with the 
HOMO and LUMO included in the active space) results are 
compared. The inset shows the energy gain of the three latter 
methods vs. RHF. In b) restricted (circles), PO- (triangles) 
and BS-unrestricted BHandHLYP (squares) results are com- 
pared; addition energies E a£ jd defined in Refs . 29 ' 30 are also 
illustrated. 
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FIG. 9: (Color online) Intramolecular charge transfer between 
the donor and acceptor moieties in TCNQ-PHP-TTF induced 
by sparkle charges as obtained by a) ab initio wavefunction 
and b) hybrid DFT with a 6-31G** basis set. In a) restricted 
HF (circles), PO-unrestricted HF (triangles), and CASSCF 
(diamonds, with the HOMO and LUMO included in the active 
space) results are compared. In b) restricted (circles) and 
PO-unrestricted BHandHLYP (triangles) data are compared. 
The addition energies E add defined in Refs.^ anc-22 are also 
illustrated. In both a) and b), the BS-unrestricted initial 
guess leads to the closed-shell restricted solution. 



correct and general approach to the problem, a single- 
determinant spin-unrestricted Hartree-Fock (UHF) ap- 
proach can provide a fair description, unlike RHF which 
produces qualitatively wrong results^. A spin-polarized 
solution allows initially paired electrons to sit in differ- 



ent spatial orbitals upon dissociation; however, this ad- 
ditional degree of freedom comes at the cost of a break- 
ing in the symmetry of the many-electron wavefunction 
resulting in spin contamination, which is sometimes re- 
ferred to as the symmetry dilemma in the literature 52 . 
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Spin-polarized (both pure and hybrid) DFT schemes are 
also known to describe dissociation problems quite satis- 
factorily^^. Their performance in describing open-shell 
molecular systems, such as singlet biradicals or dissoci- 
ating molecules, classified as type-I and type-II systems, 
respectively, have also been explore d 20 ' 54 ' 55 . 

It is therefore of interest to investigate whether single- 
determinant spin-unrestricted DFT and HF methods are 
able to deliver a step-like one-electron jump between 
weakly coupled donor and acceptor moieties upon charg- 
ing by an external field, which is the crucial feature ab- 
sent in closed-shell descriptions. For this purpose, we 
focus again on the two molecules introduced in the pre- 
vious Section. 

It is well known with both UHF or spin-polarized DFT 
that, when the starting wavefunction or charge density 
in a self-consistent field (SCF) calculation is built in a 
conventional manner, the converged solution is always 
identical to the closed-shell result, even for cases where 
this solution is physically unstable. This is because the 
original symmetry of the initial guess cannot be broken 
by the self-consistency cycles, implying that if the spatial 
orbitals are initially the same for both spins, this will be 
also the case at the end. Thus, the spin-symmetry of 
the initial guess has to be broken in some way in order 
to reach the energetically lowered and correct open-shell 
solution for systems where such a solution is likely to 
exist. 

Two ways for achieving this are commonly used 2 ^: 
(i) occupied and vacant orbitals are permuted for one 
spin only (permuted orbital scheme PO) thus creating 
directly an open-shell initial guess (this is implemented 
in Gaussian with the keyword guess=alter); or (ii) the 
HOMO and LUMO of the closed-shell initial guess are 
mixed (broken symmetry scheme BS), thereby breaking 
its spatial symmetry and introducing in the new guess 
some two-configurational character (this is implemented 
in Gaussian with the keyword guess=mix). For both 
schemes, some degree of spin contamination due to the 
mixture of singlet and triplet states before the SCF cycle 
is expected to be also found in the converged open-shell 
result. Nevertheless, the PO initial guess can offer a pos- 
sibility to treat type-I systems at a single-determinant 
level while the BS initial guess can lead to reasonable 
results for type-II systems^ -. 

Since the closed-shell HOMO and LUMO of the D- 
A system in our calculations correspond mainly to the 
HOMO of the donor and the LUMO of the acceptor, 
respectively, permuting or mixing them enforces charge 
transfer between these moieties in the initial guess and 
thus favor the desired solution. We stress that with both 
methods the initial guess engineering does not lead to 
an open-shell solution different from the closed-shell one 
in the absence of an external field (Figs. [5] and [5]). For 
both molecules, DFT calculations with LDA or B3LYP 
functionals collapse into the closed-shell solutions what- 
ever the initial guess. We thus focus hereafter on results 
obtained with Hartree Fock and hybrid DFT using the 



BHandHLYP functional. 

For the molecule with partial overlap between donor 
and acceptor orbitals (TCNQ-PP-TTF in Fig.©, we find 
that both approaches to define the initial guess in open- 
shell systems (PO and BS) allow us to recover distinct 
steps in the electron transfer with UHF (Fig. [5^,) and 
spin-polarized hybrid DFT (Fig. [HJd) . There are, how- 
ever, significant differences in the stability range of the 
UHF results obtained with the two approaches compared 
to CASSCF, providing the correct benchmark solution 
with the lowest total energy (inset of Fig. [5^,). Some 
lowering of the CASSCF energies by a constant amount 
for and 2 electron transfer (regions I and III in Fig. [3]) 
with respect to RHF is due to the stabilization arising 
from closed-shell dynamic electron correlation. More sig- 
nificant is the CASSCF energy lowering in the region 
II, where the RHF solution is inadequate; the quality 
of the UHF solutions is determined by the proximity of 
their energies to the CASSCF energy curve. As can be 
seen from Fig. 8a and its inset, PO-UHF tends to devi- 
ate from CASSCF by inducing the electron transfer too 
early, which can be intuitively understood by the fact 
that, within this approach, the initial guess is already 
rather close to the solution for the full-electron trans- 
fer configuration. The BS-UHF method, with the initial 
guess configuration quite close to RHF, has the opposite 
problem so that the electron jump occurs for a larger field 
compared to CASSCF. 

For both approaches, the UHF solutions are reasonable 
approximations to CASSCF only in the narrow range 
of field where their total energies coincide (see inset of 
Fig. [8^) though they are even there spin contaminated 
with (S 2 ) ~ 1.3. The open-shell solutions from spin- 
polarized DFT with the BHandHLYP functional for the 
same molecule (Fig. [5})) do not appear to suffer from the 
same deficiencies. In this case, both PO and BS con- 
verge to similar solutions that are also quite close to the 
CASSCF benchmark in Fig. Hi. 

For the non-conjugated TCNQ-PHP-TTF molecule 
(Fig. O, the BS-UHF solution coincides with RHF 
(Fig. [5^). Since there is here no overlap between the 
HOMO and LUMO levels which are localized on the 
rather distant D and A moieties, the assumption of BS- 
UHF that only small corrections to the RHF initial guess 
would be sufficient to direct the open-shell solution sim- 
ply fails. On the other hand, PO-UHF (Fig.[5Ji) and spin- 
polarized DFT with a BHandHLYP functional (Fig. [9b) 
do recover a step-wise electron transfer. These solutions 
are not spin pure and point to singlet-triplet mixtures 
with (S 2 ) of about 1.3 and 1, for UHF and spin-polarized 
DFT, respectively. For a quantitative assessment of the 
usefulness of our results, we need to compare them to the 
correct pure singlet state described by CASSCF for the 
whole range of external field. Fig.[5Ji shows that the PO- 
UHF step occurs at a field threshold much lower than 
for CASSCF (in close similarity with Fig.[5k), whereas it 
occurs at the onset of partial electron transfer in RHF, 
as explained in Section IIIII The behavior of the spin- 
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polarized DFT results (Fig. [BJd) appears to be rather sim- 
ilar to PO-UHF (Fig. EJa) for TCNQ-PHP-TTF, which 
puts it at odds with the CASSCF reference. 

We emphasize again that our results obtained with 
LDA or BLYP as well as with the quite commonly used 
hybrid functional B3LYP (which we do not show here), 
exhibit artificial continuous charging even when allowing 
for spin-polarization. Whatever the choice of the ini- 
tial guess, with these parameterizations of the XC func- 
tional, the open-shell ansatz collapses in general wrongly 
into the closed-shell solution for the TCNQ-PP-TTF 
molecule. In the absence of any orbital overlap, which is 
the case e.g. for TCNQ-PHP-TTF, we encounter spuri- 
ous behavior in the SCF cycle and convergence problems. 
It is only with BHandHLYP (50% of HF exchange) that 
we can observe the physically correct transfer of integer 
amounts of electrons. 

As a summary for this subsection, we point out that 
spin-polarized hybrid DFT calculations both with PO 
and BS initial guesses appear to be a viable solution 
for describing stepwise electron transfer with quite good 
threshold values for the external field, at least for the 
TCNQ-PP-TTF molecule with a strong conjugation be- 
tween the D and A units. For TCNQ-PHP-TTF charac- 
terized by a smaller coupling, the transfer steps are cor- 
rectly described at a qualitative level whereas the thresh- 
old field values deviate quite substantially from the cor- 
rect reference provided by CASSCF. The performance 
of UHF are quantitatively not very convincing for both 
molecules although steps are also observed in the curves 
with the BS and PO approaches. This failure of a Hartree 
Fock ansatz is a further evidence that the problems en- 
countered with single-determinant closed-shell methods 
for the description of electron transfer in the weak cou- 
pling regime are not necessarily intimately linked to the 
self-interaction issue in DFT. 
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D. The role of the HOMO-LUMO gap in electron 
transfer processes 

The HOMO-LUMO gap in a closed-shell Hartree-Fock 
description is an important parameter that fully deter- 
mines the thresholds or onsets for electron transfer pro- 
cesses. Our analytical model introduced in Sec. IIIII pre- 
dicts, in the absence of overlap between the frontier or- 
bitals, that the Coulomb integral is equal to the energy 
difference between the two levels at the threshold value 
of the external field (see Eq. [5]). We now address here 
whether our model predictions are reflected in the evo- 
lution of the HOMO-LUMO gap with the external field, 
as evaluated from first-principle calculations both in the 
framework of HF and DFT with hybrid XC function- 
als. We show in Figs. [TOh. and b for the TCNQ-PP-TTF 
and TCNQ-PHP-TTF molecules, respectively, the evolu- 
tion of the HOMO-LUMO gap provided by closed-shell 
RHF calculations and hybrid DFT with a BHandHLYP 
functional as a function of the external field. The gap 



FIG. 10: (Color online) HOMO-LUMO gap as a function of 
the sparkle charge for a) TCNQ-PP-TTF and b) TCNQ-PHP- 
TTF at the closed-shell Hartree-Fock (open symbols) and hy- 
brid DFT BHandHLYP (closed symbols) levels. 



decreases with a growing external field up to a certain 
threshold at which electron transfer is initiated and then 
remains constant until the two electrons have been ex- 
changed between the frontier orbitals, in full consistency 
with the analytical model. In close similarity with the 
differences observed between the charging curves of the 
two molecules in Figs. [5] and [91 the shape of the curves 
in Fig. rrOh, for TCNQ-PP-TTF is found to be somewhat 
smoothened due to orbital overlap when compared to the 
strictly linear behaviour obtained for TCNQ-PHP-TTF 

in Fig. rrub. 

For the threshold gap Ae c , which we define as the size 
of the HOMO-LUMO gap at the onset of orbital mixing, 
we obtain values of 2.26 and 1.03 eV from HF and 1.03 
and 0.52 eV from hybrid DFT for TCNQPP-TTF and 
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TCNQ-PHP-TTF, respectively, from Fig. \M In Fig. H 
we showed that the asymptotic limit of this threshold 
evolves as 1/d for large d within the framework of HF. 
Making use of this linear relationship we can formally cal- 
culate from these critical gaps effective electron transfer 
distances as 

^f(m = ^ Wry (ID 

where atomic units have been assumed and related con- 
version factors have been omitted. In this way we de- 
rive d eff to be 6.4 A for TCNQ-PPTTF and 13.9 A for 
TCNQ-PHP-TTF, which roughly corresponds to the sep- 
aration between the centers of the TTF and TCNQ moi- 
eties (5.9 A for the PP and 14.2 A for the PHP bridge, 
respectively) . 

Figs. [TOk and b also show that, when moving from 
a pure HF description to DFT with a BHandHLYP XC 
functional containing 50% of HF exchange, Ae c decreases 
by the same factor, i.e., by one half for TCNQ-PHP-TTF 
and a close value for TCNQ-PP-TTF. This is not obvious 
to explain since the meaning of the electronic eigenener- 
gies of MOs is different within HF and DFT. In pure 
DFT (with semilocal exchange) these eigenvalues corre- 
spond to chemical potentials and charge transfer should 
only take place when the energies of the frontier orbitals 
are actually crossing due to the shift induced by the ex- 
ternal field. In HF, on the other hand, charge transfer 
sets in when they are still separated by a threshold gap 
compensating the Coulomb integral between the two or- 
bitals as we have shown in Section UTIK . 

The virtually linear relation between the threshold MO 
gap and the percentage of HF/DFT exchange obtained 
from Fig.fTOlnow allows us to generalize Eq.[TT1 derived so 
far within HF by introducing the fraction of HF exchange 
in the XC functional of hybrid DFT Jhf-x as a variable. 
We then obtain: 



d eff 



E. Practical schemes for calculating addition 
energies from closed-shell DFT calculations 

We have demonstrated above that the qualitatively 
correct step-like behavior of charging curves can be ob- 
tained from the open-shell (spin-polarized) solutions of 
DFT with hybrid functionals, though the quantitative 
values for the onsets still differ from the CASSCF refer- 
ence. There is, however, a different way to approach the 
problem within DFT, where it becomes possible to de- 
rive quantitatively realistic values for CB addition ener- 
gies from closed-shell calculations although the charging 
curves are qualitatively wrong. 

In order to explain this apparent contradiction it has 
to be emphasized that the HOMO and LUMO in a single- 
particle KS scheme does not match in general the total- 
energy difference between the ground state and lowest 
charged states when the size of the HOMO-LUMO gap 
is finite 59 . This has been recently exploited for realis- 
tic calculations of Fi a dd with standard DFT techniques in 
three different ways: i) for metal particles of finite size, 
a modified KS gap has been suggested, where the ener- 
getic difference of the LUMO for charged and uncharged 
clusters has been directly taken into account^; ii) for 
the description of the gap at Cgo-metal interfaces, the 
charging energy has been obtained by using a constrained 
DFT formalism® 1 , where the occupation of hand-picked 
orbitals can be defined as a constraint in the input 62 ; iii) 
within a NEGF-DFT framework, E a( M has been defined 
via threshold values of an external gate voltage V gate de- 
termined via a midpoint integration rule from induced 
charge transfer between small molecules (H2 and ben- 
zene) and lithium wires 29 ; this method has been also ex- 
tended to aluminum surfaces and was shown to correctly 
describe screening effects^. According to the midpoint 
integration rule in Rcfs ; 29 i 30 , the onsets can be deter- 
mined from the voltages required to move 0.5 and 1.5 
electrons, respectively. As a matter of fact, Figs. |S] and [S] 
show that the onsets determined from the closed-shell 
DFT solutions match almost exactly the CASSCF re- 
sults. 



which is valid for HF (/#f-x=1), DFT with semilo- 
cal exchange (/hf-A'=0) and by design for hybrid DFT, 
where Jhf-x is defined by the choice of XC functional. 

For all discussed methods, the threshold external po- 
tential for electron transfer Ae c is physically defined by: 
i) the energy gap between the relevant donor and accep- 
tor orbitals of the system without the perturbation of the 
external field; and ii) the slope of the MO energies as a 
function of the field which is purely electrostatic in its 
origin 5 ^. Since the HOMO-LUMO gap (or bandgap in 
solids) is severely underestimated by DFT with local or 
semi-local XC functionals but, on the other hand, over- 
estimated by HF, hybrid DFT schemes provide in many 
cases a solution which better agrees with experimental 
data££. 



V. SUMMARY 

In conclusion, we have presented a configuration in- 
teraction description of electron transfer between weakly 
coupled organic molecules. The use of a multidetermi- 
nant wavefunction approach points to the key role of 
many-body effects in such a scenario. Our approach 
yields the physically correct step-like features while a 
closed-shell ansatz, either in the DFT or RHF framework, 
introduces severe limitations that we have explained via 
an analytical derivation. First-principle calculations in 
the HF and DFT framework corroborate our predictions 
for the deficiencies of closed-shell solutions. In a proper 
open-shell treatment, step-like jumps that are the hall- 
mark of charge quantization can be recovered with UHF 
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and spin-polarized DFT with a hybrid exchange. Fur- 
thermore, by relating the onset of charge transfer to the 
zero bias HOMO-LUMO gap, we discussed its origin and 
meaning in the context of both HF and DFT. Our work 
has also been connected to recently proposed practical 
schemes for calculating the addition energies in electron 
transport experiments in the CB regime from closed-shell 
DFT calculations. 
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